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(Nl Abstract 

This paper revisits an algorithm for isolating real roots of univariate poly- 
C/3 nomials based on continued fractions. It follows the work of Vincent, Uspen- 

sky, Collins and Akritas, Johnson and Krandick. We use some tricks, espe- 
cially a new algorithm for computing an upper bound of positive roots. In 
i — i this way, the algorithm of isolating real roots is improved. The complexity of 

our method for computing an upper bound of positive roots is 0(n log(w+l)) 
where u is the optimal upper bound satisfying Theorem [3] and n is the de- 
q gree of the polynomial. Our method has been implemented as a software 

package logcf using C++ language. For many benchmarks logcf is two or 
i— i three times faster than the function Rootlntervals of Mathematica. And 

^ it is much faster than another continued fractions based software CF, which 

seems to be one of the fastest available open software for exact real root isola- 
tion. For those benchmarks which have only real roots, logcf is much faster 
than Sleeve and eigensolve which are based on numerical computation. 

Keywords: Univariate polynomial, real root isolation, continued fractions, 
i—i computer algebra 

> 
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1. Introduction 

Real root isolation of univariate polynomials with integer coefficients is 
one of the fundamental tasks in computer algebra as well as in many appli- 
cations ranging from computational geometry to quantifier elimination. The 
problem can be stated as: given a polynomial P G Z[x], compute for each 
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of its real roots an interval with rational endpoints containing it and being 
disjoint from the intervals computed for the other roots. The methods of 
isolating real root can be divided into three kinds. The first kind consists 
of the subdivision algorithms using counting techniques based on, e.g., the 
Strum theorem or Descartes' rule of signs. This kind of methods count the 
sign changes (of Sturm sequence or coefficients of some polynomials) in the 
considered interval and if the sign changes reach 1 or 0, the procedure re- 
turns from this interval. Otherwise it subdivides the interval and compute 
recursively. The symbolic implementations of these methods can be found in 
[3J HH] and the symbolic-numberic algorithms implementations can be found 

m data Bang. 

The second kind takes use of the continued fraction (CF) algorithms 
[U E21 120] . These methods are highly efficient and competitive [T8| El [11] . 
Especially, [11] provides a test datasets consisting of 5000 polynomials from 
many different settings, with results indicating that there is no best method 
overall. However one can say that for most instances the solvers based on CF 
are among the best methods. In this paper we modify a real root isolation 
algorithm based on CF method to obtain a more efficient tool logcf . 

The third kind is based on Newton-Raphson method and interval arith- 
metic. The search space is subdivided until it contains only a single real root 
and Newton's method converges. When the polynomial is sparse and has 
very high degree, this method will be much faster than other methods. The 
symbolic implementations of this kind of methods can be found in (23J EI] 
and the numeric implementations can be found in [151 UH]- 

Those methods based on CF compute the continued fraction expansion 
of the real roots of a polynomial in order to compute isolating intervals for 
real roots. One important step is the computation of upper bounds of the 
positive real roots of some polynomials. There are several classic methods to 
compute such upper bounds, such as Cauchy bounds, Lagrange-MacLaurin 
bounds and Kioustelidis' bounds. There are many recent works about the 
upper bound of the positive roots of univariate polynomials [TJl EH El El H] • 
Some methods for computing such bounds are of 0(n) complexity but the 
results are very coarse like Cauchy bounds. Some methods are of 0(n 2 ) 
complexity but their bounds are sharper such as the method presented in 
[I]. The balance between precision and effect for computing such upper 
bounds has to be taken into account. 

We provide a new method for computing such bounds with time complex- 
ity 0(n log(w-l-l)), where u is the optimal upper bound satisfying Theorem|3} 
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Besides, compared with |4j, when Algorithm [5] return true (the upper bound 
is less than 1), our upper bound is at most two times that in In this way, 
the algorithm of isolating real roots is improved. Our method has been im- 
plemented as a software package logcf using C++ language. For many bench- 
marks logcf is two or three times faster than the function Rootlntervals of 
Mathematica. And it is much faster than another continued fractions based 
software CF, which seems to be one of the fastest available open software 
for exact real root isolation. For those benchmarks which have only real 
roots, logcf is much faster than Sleeve and eigensolve which are based 
on numerical computation. 

The rest of this paper is organized as follows. Section 2 reviews the 
main algorithm for real root isolation based on CF. Section 3 presents a 
new algorithm for computing an upper bound of positive roots. Section 4 
lists some tricks used in logcf. Section 5 lists the comparative experimental 
results of our algorithm and other software. 



2. Algorithm based on CF 

In this section, we first recall Descartes' rule of signs, which gives a bound 
on the number of positive real roots. Then the Vincent theorem, which can 
ensure the termination of algorithms based on CF, is presented. Finally, we 
review an algorithm of real root isolation based on CF. 

As usual, deg(p) denotes the degree of univariate polynomial p. The 
derivative of polynomial p with respect to the only variable is denoted by p' 
and gcd(/, g) means the greatest common divisor of polynomials / and g. 

Notation 1 (Sign variation). Let S = {ao, } be a finite sequence 

of non-zero real numbers. Define V(S), the sign variation of S, as follows. 



V(S) =0 if \S\ < 1 
V(a , • • • , a n _i,a n ) = 



V(a , . . . , a n -i) + 1 if a n _ia„ < 0; 
V(a , . . . , a„_i), otherwise. 



If some elements of S are zero, remove those zero-elements to get a new 
sequence and define V(S) to be the sign variation of this new sequence. 

Theorem 1 (Descartes' rule of signs). Suppose p = Y^i=o a i xl ^ ^ as m 
positive real roots, counted with multiplicity. Set V(p) = V(a , a 1; . . . , a n ). 
Then m < V(p), and V(p) — m is even. 
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Theorem 2 (Vincent's theorem). Let P(x) be a real polynomial of degree n 
which has only simple roots. It is possible to determine a positive quantity 
5 so that for every pair of positive real numbers a and b with \b — a\ < 5, 
the coefficients sequence of every transformed polynomial of the form P(x) = 
(l + x) n P{^) has exactly or 1 sign variation. The second case is possible 
if and only if P(x) has a single root within (a,b). 



Algorithm 1. main 



Input: A non-zero polynomial P(x) G 

Output: /, a set of real root isolating intervals of P(x). 

1 I = 0; 

2 if deg(P) == then 

3 ^ return J; 

4 F = gcd (p iP /) ; /* square free */ 

5 if P(0) == then 

6 
7 



Ladd([0,0}); /* add [0,0 



to set / */ 
dec(P); /* Algorithm [2 */ 

8 J.addAll(cf(P)); 

/* add all the positive root intervals to set / */ 
/* cf is described as Algorithm [4] */ 

9 p = -p; 

10 PaddAll(cf(P)); 



CF based procedures will continue subdividing the considered interval 
into two subintervals and make a one to one map from (a, b) to (0, +00) by 
P(x) = (1 + x) n P(^f) until V(P) equals 1 or 0. Therefore, Theorem g 
guarantees the termination of these procedures. 

Definition 1. As in [2], we define the following transformations for a uni- 
variate polynomial P(x). 

P(P(x))=x"(P(i)), 

H x (P(x)) = P(Xx), 
T(P(x)) = P(x + 1). 
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T(P) is also called Taylor shift one [HI E]. In our experiments when 
Algorithm [6] is used for computing upper bounds, T(P) takes more than 
ninety percent of running time^J We have considered methods in [9] for 
computing T(P), but finally we chose the classical method (Horner's method) 
for its simplicity. In future work we will use Divide & Conquer method 
which is the fastest in [9]. We think this substituting will still improve the 
performance of our method. 

Algorithm 2. dec 
Input: P = a n x n + a n _\X n ~ x + • • • + a\X + a ,n > 0. 
Output: P = a n x n ~ l + a n _iX™~ 2 + • • • + a,2X + ai . 



Algorithm 3. loglb 
Input: P~eZ\x\. 

Output: rootJb, a lower bound of positive roots of P. 

1 P = R(P); 

2 root Jb =logup(P); /* logup is described as Algorithm 6 */ 



{(min { — , ^1 , max {-,-}) if cd ^ 0; 
yea) yea} 
(0, oo), otherwise. 

Using the above notations and definitions, an algorithm for isolating all 
the real roots of a nonzero univariate polynomial is described as Algorithm 
[T] Algorithm |4j which has only a little modification of the algorithm in [1] , 
is presented here to make our subsequent description clearer. 

3. A new algorithm of computing upper bounds 

One key ingredient of CF based methods is the computation of upper 
bounds of the positive real roots of some polynomials. We give in Theorem [3] 
a new characteristic of such upper bounds of univariate polynomials. A new 
algorithm based on this theorem, Algorithm [6j is proposed for computing 
upper bounds of positive real roots. 



Hhe result of GNU gprof. 
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Algorithm 4. cf 



Input: A squarefree polynomial F G Z[x] \ {0}. 

Output: roots, a list of isolating intervals of positive roots of F. 

1 roots = 0; s = V(F); 

2 intstack = 0; mts£acfc.add({l, 0, 0, 1, F, s}); 

3 while intstack ^ do 
{a, 6, c, d, P, s} = intstack. pop();/* pop the first element */ 
a = loglb(P); 
if a > 1 then 

{a, c, P} = {«a, ac, H a (P)}; {b, d,P} = {a + b,c + d, T(P)}; 
if P(0) == then 
Lroat a .add([5,5]);P=f; 

s = V(P); 
if s == then 
^ continue; 

else if s == 1 then 

j roots. &dd(intvl(a, b, c, d)); continue; 

{Pi, a-i, &i, ci, di, r} = {T(P), a, a + 6, c, c + d, 0} 
if Pi(0) == then 
Lroot S .add([|,|]);P 1 = f;r = l; 

si = V^(Pi); {s 2 , a-2, &2, C2, d 2 } = {s - si - r, 6, a + 6, d, c + d}; 
if s 2 > 1 then 

P 2 = ( x + l) dc s(-P)T(P); 
if p 2 (0) == then 



P 2 



; s 2 = V(P 2 ); 



if si == 1 then 

| roots. add(intvl(ai, bi, c±, di)); 

else if si > 1 then 

| mtstac/c. add({ai, &i, ci, d±, P±, Si}); 

if s 2 == 1 then 

j roots. add(intvl(a 2 , b 2 , c 2 , d 2 )); 

else if s 2 > 1 then 

|_ mtstac/c.add({a 2 , 6 2 , c 2 , d 2 , P 2 , s 2 }); 
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Theorem 3. Suppose P = a n x n + a n _ 1 x n 1 + - ■ - + aiX + ao (a n > 0) is a uni- 
variate polynomial in x with real coefficients. Then a nonnegative number u is 
an upper bound of positive roots of P ifu satisfies min™ =0 |XX=j a i ut ~''^ — 0- 

PROOF. If n == 0, then P is a nonzero constant and any positive number 
is its upper bound of positive roots. 

Otherwise, if b > u, we claim that Y^=j a ft~^ > ^i=j a i u% ~^ f° r an J 
j = 0, ...,n- 1. 

When j = n-l, E^n-i ^ n+1 - Etn-i a t u l ~ n+1 = a n (b-u)> 0. The 
claim holds. 

Assume the claim holds when j = k. When j = k — 1, Yl^k-i a ib l ~ k+1 = 
(ZlILfc a ^ i_fc ) b + a k-i- B y assumption Yn=k a ib l ~ k > Y^=k a i ul ~ k ^ °- 
Since b > u > 0, (Etk^ k ) b > (J2i=k a i ui ~ k ) u and Etk-i a '^ k+1 > 
E"=fc-i aiU l ~ k+l - So J27=j a i bt ~ j > EILj a i ui ~ j for any j = 0, . . . , n - 1. 

By the above claim, P{b) = ^™ =0 aft > when b > u. Because b is 
arbitrarily chosen, u is an upper bound of the positive roots of P. 

The following theorem was given by Akritas et al. in [31 H], which com- 
putes positive root upper bounds of univariate polynomials. 

Theorem 4 (Akritas-Strzebohski-Vigklas, |3j). Let P(x) = a n x n +a n _ix n ~ 1 + 
■ ■ ■ + do {a n > 0) be a polynomial with real coefficients and let d(P) and t(P) 
denote the degree and the number of its terms, respectively. 
Moreover, assume that P(x) can be written as 

P(x) = qx(x) - q 2 {x) + q 3 (x) - q±(x) H h q 2m -i{x) - q2 m (%) + g(%) (1) 

where all the coefficients of polynomials qi(x) (i = 1, 2, . . . , 2m) and g(x) are 
positive. In addition, assume that for % = 1, 2, . . . , m we have 

q 2i -i{x) = c 2i -i,ix e2 >- 1 ' 1 + ■■■ + c^.i,^.^ 621 - 1 ^-! 

and 

q 2l (x) = b 2iil x e2 ^ + ■■■ + b 2lMi x e2 ^ 

where e 2i _i,i = d{q 2i -i), e 2i ,i = d(q 2i ), t 2i _ x = t(q 2i -x), and t 2i = t(q 2i ) and 
the exponent of each term in q 2 i-\(x) is greater than the exponent of each 
term in q 2 i(x). If for all indices i = 1, 2, . . . , m, we have 

t{q2i-i) > t(q 2 i), 
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then an upper bound of the values of the positive roots of p(x) is given by 




up = max { max { I — } /> /> (2) 

i=l,2,...,m 1 j=l,2,...,t 2i 11 ' 11 

for any permutation of the positive coefficients C2i~ij, j = 1, 2, . . . , £21-1- 
Otherwise, for each of the indices i for which we have 

^2i— 1 < t2i, 

we break up one of the coefficients of fe-i^) into t 2 i — £21-1 + 1 parts, so 
that now £((?2i) = t(<72i-i) and apply the same formula given above. 

We shall show in Theorem [6] that the bound given by Theorem [3] is better 
than that given by Theorem |4j 

Theorem 5. Let P(x) = a n x n + a„_ix n_1 + ■ ■ ■ + a (a n > 0) be a polyno- 
mial with real coefficients and u denote an upper bound of positive roots of p 
obtained by Theorem^ then min£ =0 {Y^=k a i ul k } > 0. 

Proof. For every Oj < 0, by Theorem [IJ there exist c^x 6 ^ and &j 2 x e<2 , 
respectively, such that > e^ and c^w 6 ^ > bi 2 u e ^. By Theorem [4] 6 i2 x e ^ is 
the term —aiX 1 and Ci x x e ^ is either a whole or a part (broken up by Theorem 
111) of a positive term of p. 



For every a, > 0, by Theorem 



(EaKO.e^-Cu) < So ElLfc ^ > 



Zir=fc,a l <o ( c ii Men - k^ 6 ' 2 ) > for an Y = 0, 1, . . . , n. Then £)" =fc a;w* fc > 
for any = 0, 1, . . . , n and min^ =0 \Y™ =h a i u t ~ k ^ > 0. 

Theorem 6. Let P(x) = a n x n + a n -\X n ~ l + ■ ■ - + a (a n > 0) be a polynomial 
with real coefficients. Let U\ denote the optimal upper bound of positive real 
roots satisfying Theorem^ and U2 denote the optimal upper bound of positive 
real roots satisfying Theorem \Q then U\ < U2 and the strict inequality can 
hold. 

Proof. By Theorem [5j m <u 2 . 

Let P(x) = x 2 + x — 2. Then w 2 — \/2 an d u% = 1. So u\ < U2 for this P. 

Theorem 7. Let P(x) = a n x n + a n -ix n ~ l H h a (V(P) > 0) 6e a 

polynomial with real coefficients. Let u denote the output of Algorithm^ and 
U\ denote the optimal upper bound of P satisfying Theorem [3| When u is 
less than or equal to 1, u < 2u\. 
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Algorithm 5. lessOne 



Input: P = a n x n + a n ^ix n ~ l + • • • + d\X + a G Z[x], 3oj, a„aj < 0. 
Output: true: the positive root bound of P must be less than 1; 

false: cannot determine whether the bound is less than 1. 

1 start = n — 1; 

2 lastNeg = 0; 

3 hSign = sign(a n ); 

4 while sign(a/ astA r efl ) + hSign ^ do 

5 j lastNeg = lastNeg + 1; 

6 while sign(a stari ) + hSign ^ do 

7 | start = start — 1; 

8 cfSum = abs(a„); 

9 i — n — 1; 
10 j = start; 

n /ast = start; 

12 while i > lastNeg — 1 and j > lastNeg — 1 do 

13 
14 
15 

16 
17 

18 
19 

20 
21 
22 

23 
24 

25 
26 



if sign(c/S'-um) < then 

while i > last and sign(aj) ^ hSign do 
\_i = i-l\ 

if % == last then 
j return false; 

cfSum = cfSum + abs(aj); 
i — i — 1; 

else 

if j == lastNeg — 1 then 
] return true; 

while j > lastNeg and sign(aj) + hSign ^ do 

i - .7 - 1 : 
cfSum = cfSum — abs(aj); 
last = j; 



27 return true; 
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Algorithm 6. logup 

Input: P = a n x n + a n -\x n ~ x + • • • + a\x + ao G Z[x], 3aj, a n aj < 0. 
Output: an upper bound of the positive roots of P. 

1 start = n — l; lastNeg = 0; hSign = sign(a n ); 6ase = 1; 

2 if -dessOnefP,) then 

3 i return 2; 

4 while sign(ai ast Ne 9 ) + /iSigra 7^ do 

5 I lastNeg = lastNeg + 1; 

6 while sign(a sta rt) + hSign ^ do 

7 j start = start — 1; 

8 i = n; 

9 while % == n do 

i — n — 1] 
j = start; 
cfSum = abs(a n ); 

while i > lastNeg — 1 and j > lastNeg — 1 do 
if sign(c / Sum) < then 

while i > j and sign(aj) 7^ hSign do 
j % = % - 1; 

if % == j then 

|_ break; 

cfSum = cfSum + ahs(ai)2^ base ; 
i — i — 1; 

else 

if j == lastNeg — 1 then 
j j = lastNeg — 2; break; 

while j > lastNeg and sign(aj) + hSign 7^ do 
J-J-l: 

c/S«m = c/S«m - abs(ai)2^-^ base ; 

3=3- 1; 



10 
11 
12 

13 
14 
15 
16 

17 
18 

19 
20 

21 
22 
23 

24 
25 



26 
27 



28 
29 



30 return 



if j == lastNeg — 2 then 
j base = base + 1; i — n\ 

1 
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Proof. In Algorithm 



if 2^ > M i> tnen min i=o {Ta=j a i {ztk^y J } > 
by the proof of Theorem 3] and thus the loop does not terminate at this step. 
So when Algorithm [6 returns, base must satisfy < u\. Therefore, the 
output u = 2 ,, as 1 c _i and u < 2u±. Obviously, this algorithm will terminate. 

Furthermore, min"~Q |Xir=j a « ( 2^-1 T ^} — ^y Theorem jsj So, u = 
26a s e _i is an upper bound of p. 

Corollary 8. Let P(x) = a n x n + a n „ x x n - x + • • • + a (V(P) > 0) be a 
polynomial with real coefficients. Set u to be the optimal upper bound of 
positive roots of P satisfying Theorem [3| Then Algorithm [6| costs at most 
0{n\og{u + 1)) additions and multiplications. 

4. Tricks 

Variable substitution If P(x) £ Z[x] and P(x) = P\(x k ) (k > 1), then 
substitute y = x k in P. Obviously, deg(Pi,y) = deg ^ P '^ . We first isolate 
the real roots of Pi then obtain the real roots of P. We can see in Figure 
[2] that degree is a key fact affecting the running time. Using this trick, we 
can greatly reduce the running time of ChebyshevT and Chebyshevll when 
each term of the polynomials is of even degree. The running time on such 
polynomials can be found in Table [2] The same trick was also taken into 
account in [T5] . 

Incomplete termination check If P(x) £ Z[x] and V(P) = 2, we may 
try to check whether the sign of P(l) is the same as the sign of the leading 
coefficient of P. If they are not the same, then P has one positive root in 
(0, 1) and the other one in (1, +00). So, we can terminate this subtree. Since 
the whole logcf procedure is a tree and logcf spends more than 90 percent 
of the total time on computing T(P), this trick may improve the efficiency 
of the algorithm greatly. 

5. Experiments 

5.1. Implementation 

The main algorithm for isolating real roots based on our improvements 
has been implemented as a C++ program, logcfj^J Compilation was done 



The program can be downloaded through http://www.is.pku.edu.cn/~dlyun/ 
logcf 
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using g++ version 4.6.3 with optimization flags -02. We use Singular [TD] 
to read polynomials from files or standard input and to eliminate multi- 
factors of polynomials. We use the GMFj^] (version 5.05), arbitrary-length 
integers libraries, to deal with big integer computation. All the benchmarks 
listed were computed on a 64-bit Intel(R) Core(TM) i5 CPU 650 @ 3.20GHz 
with 4GB RAM memory and Ubuntu 12.04 GNU/Linux. 

5.2. Benchmarks 

5.2.1. W n 

Wilkinson polynomials: W n = U.^ =1 (x — i). The integers 1,2, ... ,n are all 
the real roots of W n . 

5.2.2. mW n 

Modified Wilkinson polynomials: mW n = W n — 1. 

If n > 10, mW n has n simple real roots but most of them are irrational. 

5.2.3. IW n 

The distance between W n 's two nearest real roots is 1 and the distance 
between mW^s two nearest real roots is nearly 1. We construct new poly- 
nomials IW n = Hf =1 (ix — 1), which have a completely different distance 
between any two nearest real roots. 

5.2.4. mIW n 

We modify IW n into mIW n = IW n — 1 for the same purpose as we 
construct mW n . Most real roots of mIM n become irrational. 

5.2.5. T n 

ChebyshevT polynomials: T = 1, T\ — x,T n+1 = 2xT n — T n _ 1 . T n has n 
simple real roots. 

5.2.6. U n 

ChebyshevU polynomials: U — 1, U 1 — 2x, U n+ i = 2xll n — U n _\. U n has 
n simple real roots. 

5.2.7. L n 

Laguerre polynomials: L = \,L\ = l—x,L n+ i(x) = ( 2w+1 ~ x )-^W~ ,iL "-:lO e ) _ 
Obviously, n\L n is a polynomial with integer coefficients. 



3 http://gmplib.org/ 
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5.2.8. M n 

Mignotte polynomials: x n — 2 (5s — l) 2 . If n is odd, M n has three simple 
real roots. If n is even, it has four simple real roots. 

5.2.9. R{n,b,r) 

Randomly generated polynomials: R(n, b, r)=a n x n + • • • + a^x + a with 
| Oj| < b, Pr[cii > 0] = 2 and Pr[a« 7^ 0] = 1 — r, where Pr means probability. 

5.3. Results 

The root isolation timings in Tables [TJ [2] and [3] are in seconds. Most of 
the benchmarks we chose have large degrees and the timings show that our 
tool is very efficient. As a built-in Mathematica symbol, Rootlntervals 
is compared with our tool logcf . The Mathematica we use has a version 
number 8.0.4.0. For almost all benchmarks, our software logcf can be two or 
three times faster than Rootlntervals. The comparative data can be found 
in Table [TJ Table [2] and Figure [2] We also consider open software, such as CF 
|llj . which seems to be one of the fastest open software available for exact 
real root isolation. Many experiments about state of the art open software 
for isolating real roots have been done in [llj, which indicate that CF is the 
fastest in many cases. In our experiments, logcf is much faster than CF. 
The comparative result can be found in Table [3j We also compare logcf 
with numerical methods eigensolve [8J and Sleeve [TTJ. As eigensolve 
computes all the complex roots, we choose W n , mW n and IW n as benchmarks 
with degrees ranging from 10 to 90, which have only real roots. Sleeve 
computes only real roots but it has weak stability. Its output on W 30 only 
has eight real roots, which is obviously wrong. Sleeve's running time^jon W 10 
is 0.022 seconds and 0.024 seconds on W^o- m these two cases our software is 
about 7 times faster than Sleeve. We compare logcf with eigensolve and 
the results are shown in Figure [TJ At the be ginning when degree is 10, the 
time costs of logcf and eigensolve are almost equal. As degree becoming 
larger, the growth rate of our tool's consuming-time is much less than that 
of eigensolve. When degree reaches 90, logcf is about 20 times faster than 
eigensolve. 



4 When running time is very short we run every case for more than ten times and 
compute the mean. 
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10.6 


4.34 


mIW 50Q 


2.73 


1.01 


mW 10 oo 


140.9 


65.62 


mlWxooo 


32.9 


15.56 



Table 1: compare with Mathematica(l) 



Benchmark 


Rootlntervals 


logcf 


Benchmark 


Rootlntervals 


logcf 


-?100 


0.056 


0.01 


^100 


0.072 


0.02 


^200 


0.39 


0.03 


-^200 


0.60 


0.16 


^300 


1.29 


0.10 


-^300 


2.2 


0.69 


^400 


3.39 


0.22 


-^400 


5.64 


1.91 


Tboo 


7.26 


0.45 


-^500 


12.24 


4.59 


^1000 


90.8 


4.96 


-^1000 


150 


72.3 


t^ioo 


0.048 


0.01 


^2000 


1.22 


0.19 


£^200 


0.35 


0.03 


^2001 


1.22 


0.20 


£^300 


1.31 


0.09 


-^4000 


8.02 


1.79 


£^400 


3.35 


0.21 


M40OI 


7.98 


1.99 


^500 


6.95 


0.44 


^6000 


33.4 


7.73 


t^iooo 


87.5 


4.81 


^6001 


33.7 


7.82 



Table 2: compare with Mathematical) 



For randomly generated polynomials, we consider different settings of 
(n, b, r) as shown in Figure [2] For each setting (n, b, r), we generate randomly 
five instances and compute the mean of five running times. In almost every 
randomly generated benchmark our logcf is two or three times faster than 
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Rootlntervals. And We can also find that degree is the main factor affecting 
the running time. 



Benchmark 


CF 


logcf 


Benchmark 


CF 


logcf 


Wioo 


0.054 


0.01 


Woo 


0.056 


0.01 


^200 


0.23 


0.015 


Woo 


0.20 


0.015 


mW 100 


0.054 


0.025 


mIW W0 


0.14 


0.01 


mW 20 o 


40.5 


0.16 


mIW 2 oo 


2.7 


0.04 


Tioo 


0.52 


0.01 


L100 


0.80 


0.02 


T200 


4.32 


0.13 


L200 


7.50 


0.16 


U100 


0.52 


0.01 




43.52 


0.03 


U200 


4.15 


0.12 


^1200 


88 


0.05 



Table 3: compare with CF 



Compare with numerical software(W n benchmarks) 



0.2 
J 0.1 


logcf 
eignsolve 









e l r 1 — « 
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degree 

Compare with numerical software(mW n benchmarks) 
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Figure 1: compare with numerical software eigensolve 
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Compare with Mathematica (R benchmarks with b= 17951 and r=1) 



logcf 

Mathematica 


00 


500 1000 15 

degree 

Compare with Mathematica (R benchmarks with degree=1000 and r=1) 










1 2 3 4 5 6 " 

b(lo") 

Compare with Mathematica (R benchmarks with degree=1000 and b= 17951 1 
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Figure 2: R(n, b, r) benchmarks with differ setting 

from ISCAS. The authors would like to thank Steven Fortune who sent us the 
source code of eigensolve and Elias P. Tsigaridas who helped us compile 
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